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ABSTRACT 

The South China Sea (SCS), the largest marginal sea in the West Pacific Ocean, is 
separated from adjacent oceans by a chain of islands. The deepest water is confined to a 
bowel-type trench, and the maximum depth is approximately 5,000 m. Most of the existing 
studies on the seasonal and interannual variability have been based only on surface 
temperature data. However a primary need is an understanding of the three-dimensional 
thermohaline circulation. The minimum curvature with spline method was used to establish 
a three-dimensional monthly-varying gridded data from the Navy’s Master Oceanographic 
Observation Data Set (approximate 189,000 profiles), covering the area of 5°N - 25°N and 
105°E - 125°E and from the surface to 400 m depth. For temperature, profiles were binned 
into 204 monthly data sets from 1968 to 1984 (17.years). Because of the paucity of salinity 
data, salinity profiles were binned into 12 climatological monthly data sets, and the monthly 
climatological mean was computed. After the gridded data set had _ established, both 
composite analysis and the Empirical Orthogonal Function (EOF) analysis were used to 
identify the major thermohaline features. The first EOF mode accounts for 26.7% of the 
variance and represents the seasonal variation. The second EOF mode accounts for 17.7% 
of the variance and represents the interannual SCS warming/cooling phases. Furthermore, 
the P-vector method was used to invert three-dimensional velocity fields from the analyzed 
temperature and salinity data. Important dynamical processes, including the Kuroshio 
intrusion, the western boundary current (counter-current), the cross basin current (under . 


counter-current), the mesoscale eddies, and the basin gyre are identified. 





TABLE OF CONTENTS 


TORN IAT OF) COIN inno es ek ee ieee 
Po SOU erINA SEA OCEANOGRAPHY no). oe. 2. oa ae 
A. WATER MASSES ae ee. See eee | 
B. MONSOON EFFECT ON CURRENT SYSTEM ............... 
1 SUMIMCMIVIONSOONM t= — "27... .. 0... 
2 WATteIVIOHSOOL weer es ry on eee 
C, MESOSCARE EDDIES: <:: 2p Gn t.ho eh ake 
D. AN AIR-SEA FEEDBACK SCENARIO .............0....... 


Il. MASTER OCEANOGRAPHIC OBSERVATION DATA SET (MOODS) . 


IV. ESTABLISHMENT OF GRIDDED DATA ............ oes 2 anaes eae 
VaR COMPOSIEE ANALYSIS”... tee ne oe sc ccs aa 
A. ANNUAL MEAN . fe ie - 

B. SEASONAL VAREMIMON — Renee ee are re ay ee 

] RCI CRACURC jac kc nee ae, ae oe re een, 

Z SANBICY oO ee ee = 3 
VI. EMPIRICAL ORTHOGONAL FUNCTION (EOF) ANALYSIS .......... 

A THERMAL VARIABLIL IY ae... .. oe. eee ee 

B. PRINCIPAL EOPIVIODES tect seo tans pt ee ee ee 

1, TheEOF1 Mode ..........0.....0eecee essence. : 


vii 


ye The second EOF Mode ........................000000 67 


CG TEMPORAL VARIABILITIES .....2..23 = ee eee 68 

te First Principal Component oxen 69 

Z. second Prancipal Component 23 2eg eee 70 

Vil. THE SOUTH CHIN ASE A CIRCULATION (eee iG 


A. GEOSTROPHIC VELOCITY INVERTED FROM DENSITY FIELD .. 77 


B. SEASONAL VARIATION OF SOUTH CHINA CIRCULATION .... 80 


I General Desi ton 7 pater eer ee. 80 

2: Surface Cine ul atl Ocean Le Seas ye ae 80 

a Subs SUI ACCs(/ Salil eee eae ns ote 82 

4. Intermediate level (400m).. BR aio ics ee 82 

3S Kuroshio Intrusion Through the Luzon Strait .............. 83 

VULI. CONGIUSIONS Seta ee ee: 97 

A. CLIMATOLOGICAL MEAN THERMOHALINE FEATURES ...... 97 

B SEASONAL THERMOHALINE VARIABHULTY ................. 98 

C KUROSHIO INTRUSION ...... pcs) SUSIE o's a eae, nate 99 

D. INTERANNUAL SCS WARMING/COOLING PHASES ........... 99 

E MULTIEDDY STRUCTURE eee a2 100 

APPENDIX A. MONTHLY TEMPERATURE ANOMALY .................. 101 
APPENDIX B. MONTHLY SALINITY ANOMALY ........................ 161 © 

APPENDIX C. EOF 3~6 FOR TEMPERATURE ANOMALY ............ =. 2 


Vill 


APPENDIX OD, THE CIRCULATION PAIRINGS ir. oie ees ee ee 


LIST OF REFERENCES 


* © © © © # © © © «© © © © © © © © © © &# © © © © © © © © © @© @© @© © © © © © @ #@ #© © © # @ @ 


ANITA DISTRIBUTION |... ec case canes eee 





LIST OF FIGURES 


Figure |. Bathymetryofthe Southi@hinaSeaw ). sees See ee. 4 
Figure 1.2 Observational Surface Currents of the Southeast Asia Watersin April...... 5 


Figure 2.1 Distribution of water masses at the sea surface and bottom layer in the northern 
plmmenind SC4 ................. ee eee... se 14 


Figure 2.2 Mean Atmospheric Surface Streamline Analysis for summer monsoon..... 15 
Figure 2.3 Observational Summer Surface Circulation of the Southeast Asia Waters .. 16 
Figure 2.4 Mean Atmospheric Surface Streamline Analysis for winter monsoon...... 17 
Figure 2.5 Observational Winter Surface Circulation of the Southeast Asia Waters ... 18 


Figure 2.6 An Air- Sea Feedback Scenario for the SCS Warm-Core Eddy Life Cycle .. 19 


Figure 3.1 The Spatial distribution of MOODS Lemipeiaimne prenics ... 9. eee 23 
Figure 3.2 The Spatial distribution of MOODS salinity profiles .................. 24 
Figure 3.4 The temporal distribution of MOODS .......................0.-5. a5 
Figure 3.4 The temporal distribution of MOODS in 1968 ....................... oS 
Figure 4.1a Cross-sections through surfaces produced with splines in tension oe 31 


Figure 4.1b Cross-sections through continental shelf and slope produced with splines in 


SeLGHSION . . .... .ggeeeeeei ee es ee 32 
Figure 5.1a The annual mean temperature field from surface to 100 meters ......... 40 
Figure 5.1b The annual mean temperature field from 125 to 400 meters ............ 4] 


Figure 5.2a The zonal cross-section view of annual mean tempefature from Latitude 14°N to _ 
J ehe Tae ou oS, 0 cr rr rca aA care oe rein MA ee i ti 42 


XI 


Figure 5.2b The zonal cross-section view of annual mean temperature from Latitude 20°N 


tombatitide 25° Nic 2. = pac 5 ee ea 43 
Figure 5.3a The annual mean salinity field from surface to 100 meters ............. “la 
Figure 5.3b The annual mean salinity field from 125 to 400 meters................ 45 


Figure 5.4a The zonal cross-section view of annual mean salinity from Latitude 14°N to 
Latitude 19°N’ — :G kOe 0... ee... . ee. ee re 46 


Figure 5.4b The zonal cross-section view of annual mean salinity from Latitude 20°N to 
De NRE UE 2a ae A a ea ae 47 


Figure 5.5a Monthly mean temperature eae at surface from January to June... . 48 
Figure 5.5b Monthly mean temperature anomalies at surface from July to December .. 49 
Figure 5.6a Monthly mean temperature anomalies at 75 meters from January to June .. 50 
Figure 5.6b Monthly mean temperature anomalies at 75 meters from July to December 51 


Figure 5.7a The zonal cross-section view of monthly mean temperature anomalies at 19°N 
mom January to June 3.44.45. en ee eee ee, ee 32 


Figure 5.7b The zonal cross-section view of monthly mean temperature anomalies at 19°N 
HOG@eTTTCIRIV APOE COMO CT: fey tetas aaa 53 


Figure 5.8a Monthly mean temperature anomalies at 200 meters from January to June . 54 


Figure 5.8b Monthly mean temperature anomalies at 200 meters from July to December 55 


Figure 5.9a Monthly mean salinity anomalies at surface from January to June........ 56 
Figure 5.9b Monthly mean salinity anomalies at surface from July to December ...... oe 
Figure 5.10a Monthly mean salinity anomalies at 75 meters from January to June... .. 58 


Figure 5.10b Monthly mean salinity anomalies at 75 meters from July to December ... 59 


Figure 5.11a The zonal cross-section view of monthly mean salinity anomalies at 19°N from 
January to: Junes 3.0525 a, Meee eee ee ee emer ens, = om 60 


Xl 


Figure 5.11b The zonal cross-section view of monthly mean salinity anomalies at 19°N from 
UU AEGPMCCEMUCES <... 3 ¢../, 0, OE TRIE... . . Jnr, . .ccthdy . 0 eee 6] 


Figure 5.12a Monthly mean salinity anomalies at 200 meters from January to June... . 62 


Figure 5.12b Monthly mean salinity anomalies at 200 meters from July to December .. 63 


Figure 6.1 EOF] Temperature anomalies ...................0 00-00 c ee eee. 74 
Bigute ©.2 EOFZ Temperature anomaliesis sss se ss. 2 76 
Figure 6.3 The first six principal components for 1968-1984 .................... 78 
Figure 7.1 Thermal wind relation at two different depth and B spiral turning angle ... . 87 
Figure 7.2a Surface Circulation from January to June.......................... 88 
Figure 7.2b Surface Circulation form July and December ....................... 89 
Figure 7.3a Sub-Surface (75 meters) Circulation from January to June............. 90 
Figure 7.3b Sub-Surface (75 meters) Circulation form July ata eccinvcte = 9] 
Figure 7.4a Intermediate Level (400 meters) Circulation from January to June ... ee 
Figure 7.4b Intermediate Level (400 meters) Shaikh form July and December .... 93 


Figure 7.5 Temperature and Salinity data from surface to 400 meters used for the 


computation of east-westgeostrophic velocity along 120.5°E .................... 94 
Figure 7.6a The monthly zonal velocity in 120°E from January to June ............ 95 
Figure 7.6b The monthly zonal velocity in 120°E form July and December.......... 96 
Figure 7.7 The Mass Transport (Sv) at the Luzon strait ....................... 97 


Xill 





ACKNOWLEDGEMENTS 


I would like to acknowledge the unfailing support of C. W. Fan and Yuchun Chen for 
their programming assistance. Special thanks to Dr. Peter C. Chu for his guidance and 
patience and to Dr. Pierre-Marie Poulain for his valuable comments. This work was funded 
by the Naval Oceanographic Office, the Office of Naval Research NOMP Program, and the 


Naval Postgraduate School. 








I. INTRODUCTION 


The South China Sea (SCS), the largest marginal sea in the West Pacific Ocean, is 
separated from adjacent oceans by a chain of islands. It contains the broad shallow Sunda 
Shelf on the south, a shelf extending from the Gulf of Tonkin to the Taiwan Strait on the 
north, a deep sea basin in the center, and the Dangerous Ground with numerous reef islands 
over the southeast (Figure 1.1). It has a bottom topography that makes it a unique semi- 
enclosed ocean basin that is overlaid by a pronounced monsoon surface wind. Extensive 
continental shelves (less than 100 m deep) are found on the western and southern parts, while 
steep slopes with almost no shelves are found in the eastern part of SCS. The deepest water 
is confined to a bowel-type trench and the maximum depth is approximately 5,000 m. 

The SCS's physical oceanographic condition has been investigated for several decades. 
Studies show that eddies in the SCS are predominantly cyclonic in winter and anti-cyclonic 
in summer, with sizes ranging from small to meso-scale, and are both cold and warm [Huang 
et al., 1994]. A survey by Wyrtki [1961] revealed complex temporal and spatial features of 
the surface currents in both the SCS and the surrounding waters. By April, the northeast 
Ponsoon has decayed and the main SCS current consists of two large cyclonic eddies in the 
western SCS, one in the north and the other in the south (Figure 1.2). With the beginning of 
the southwest monsoon in May and June, a wind drift current is formed off the coast of 
Vietnam. The westward intensification of this current is clearly visible. The wide uniform - 


drift in the northern SCS shows a deflection of the current to the right of the wind. A large 


part of the water passes south of Taiwan and into the root of the Kuroshio. The southern SCS 
is occupied by an anticyclonic eddy. 

Based on limited data sets, studies show that eddies in the SCS are predominantly 
cyclonic in winter and spring and anti-cyclonic in summer, with sizes ranging from small to 
meso-scale [Huang et al., 1994]. Both cold and warm eddies exist in the SCS. Dale [1956] 
and Uda and Nakao [1972] reported a cold eddy off the central Vietnamese Coast in summer, 
and Nitan [1970] found a cold eddy located to the northwest of Luzon. Reports from the 
South China Sea Institute of Oceanology [SCSIO, 1985] indicate that in the coved SCS a 
warm-core eddy appears in summer and winter, but it is more closer to Vietnam in summer 
at the surface. In addition, a cold-core eddy was detected in the central SCS from the 29th 
December 1993 to 5th January 1994 according to analysis of the TOPEX/ POSEIDON data 
[Soong, et aii 1995]. 

sii Chu et al. [1997a] and Chu and Chang [1997] identified a central SCS 
surface warm-core eddy in mid-May from a more complete historical data set - the U.S. 
Navy's Master Observational Oceanographic Data Set (MOODS). From the composite 
analysis of the U.S. National Centers ee Environmental Prediction (NCEP) monthly SST 
fields (1982-1994), Chu et al. [1997b] found that during the spring-to-summer monsoon 
transition (March to May) a warm anomaly (greater than 1.8°C) is formed in the central SCS 
at 112°-119°30E, 15°-19° 30'N. This warm eddy appears in the central SCS (west of Luzon 
Island) during the late spring season and strengthens until the onset of the summer monsoon 


(mid-May) and then weakens and disappears at the end of May. Although its size and 


intensity varies, the SCS warm pool releases large moisture and heat fluxes into the 
atmosphere and in turn affects the monsoon circulation. 

Most of the existing studies on the seasonal eddy variability were based on the surface 
data. What is the three-dimensional thermal structure and circulation pattern of the SCS 
eddies? In this thesis, we will use the minimum curvature spline method [Wessel and Smith, 
1992] to establish a three-dimensional monthly varying gridded data from the MOODS, 
covering the area of 5°N - 25°N and 105°E - 125°E and from the surface to 400 m depth. 
For temperature, we binned the profiles into 204 monthly data sets from 1968 to 1984 (17 
years). For salinity, a monthly climatological analysis was applied because of the paucity of 
data. After the gridded data set had been established, we utilized both composite analysis and 
the Empirical Orthogonal Function (EOF) analysis identifing the major thermal features of the 
Soe eddies. Furthermore, three-dimensional velocity fields were inverted from analysis 
of the temperature and salinity data. 

The outline of this thesis is as follows: Chapter II depicts the South China Sea 
oceanography including the water masses, monsoon effect and mesoscale eddies. Chapter II] 
describes the U.S. Navy’s Master Oceanographic Observation Data Set (MOODS). Chapter 
IV describes the establishment of gridded data, and Chapter V depicts the composite analysis 
to obtain temperature and salinity seasonal variations. Chapter VI describes the EOF method 
to obtain nonseasonal variabilities. Chapter VII presents the SCS circulation using the P 
vector inverse method, major current systems and seasonal variation, and Chapter VIII 


includes our conclusions. 
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Figure 1.1 Bathymetry of the South China Sea 
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Figure 1.2 Observational Surface Currents of the Southeast Asia Waters in April 


(From Wyrtki, 1961) 





ll. THE SOUTH CHINA SEA OCEANOGRAPHY 


AS WATER MASSES 
On the western side of the Luzon Strait, the SCS water masses can be divided into 
those four types [Wyrtki, 1961]: (1) the warm, less salinity Surface Layer Water, (2) the high 
salinity (34.8-35.2 psu) Subtropical Lower Water, (3) the low salinity (34.1-34.5 psu) 
Intermediate Water, and (4) the onginal Pacific Deep Water ( 34.65-34.7 psu). The salinity 
maximum of the Subtropical Lower Water 1s approximately 150 meters deep and the salinity 
minimum of the Intermediate Water is approximately 400 to 500 meters deep. At the Luzon 
Strait these two water masses move in opposite directions and the depth of the inversion is 
about 200 meter. The center of this inversion lies relatively uniformly at a temperature 16°C. 
Wyrtki [1961] also found that during the winter monsoon season surface water 
masses from the Pacific Ocean are transported into the SCS through the Luzon Strait, chiefly 
along the western side of the.SCS in depths between 400 and 900 meters. During the summer 
monsoon season these conditions are reversed. 
Waters from the Western Pacific enter the northern SCS through the Luzon Strait and 
| after mixing and modification form distinct water masses. Continental fresh water runoff is 
also very significant. Su and Weng [1994] classified eight water masses in the northern SCS. 


They grouped these water masses into three classes: (1) runoff diluted type of Near shore 


Diluted Water Mass (F); (11) the shallow sea modified type composed of Coastal Mixed Water ~ 


Mass (4), Warm Surface Water Mass (WS), Surface Water Mass (8) and Surface-Subsurface 


Mixed Water Mass (SU); and (111) the deep sea type, composed of Subsurface Water Mass 
(U), Subsurface-Intermediate Mixed Water Mass (U/) and Intermediate Water Mass (J). 
Distribution of these water masses at the surface are nearly parallel to the coastline (Figure © 
Zale 

Both water masses WS and S originate from West Pacific Surface Water Mass. 
Various meteorological conditions modify this water mass and ultimately generate the high 
temperature and low salinity WS water mass and the low temperature, low salinity S water 
mass. Mixing modifies the West Pacific Subtropical Subsurface Water Mass to form the S u 
water mass, which is similar to the S water mass. The U water mass, distributed at 
approximately 150 m depth, retains the high salinity layer of the West Pacific water. The U/ 
water mass is the interface between the deep high salinity layer and low salinity surface layer. 
The J water mass is an oceanic type water mass with the least modification from its source 
water mass and characterizes low salinity, distributed down to 700 m. The F water mass is 
a mixture of runoff and Near shore water mass. The M water mass undergoes the largest 
modification, mixing with SU, WS and S water masses, respectively, in varying seasons. Table 


2.1 lists the characteristics of the water masses. 














[WaterMass | Feb. | = May _| __ Aug. | __Nov._ 
F TCO 25.5-27.8 29.0-30.5 
Salinity <30.0 <30.0 
M TCC) | 13.2-16.0 22.5-27.6 28.9-30.3 21.3-24.4 
Salinity 31.20-33.60 | 31.16-33.80 | 33.03-32.91 | 31.32-34.02 
WS TC) | 22.3-24.4 26.2-30.0 28.0-30.1 25.5-26.1 
Salinity | 34.08-34.67 | 33.90-34.69 | 33.03-34.40 | 33.68-33.93 
S TCC | 20.4-22.3 2PPI6 2 23.5-28.0 22.5-25.5 
Salinity 34.43-34.80 | 33.85-34.62 | 33.65-34.62 | 33.97-34.60 
SU TCC) | 17.5-20.4 19.8-22.2 21.0-23.5 19.0-22.5 
Salinity 34.18-34.75 | 34.16-34.71 | 34.32-34.70 | 34.42-34.76 
U TCO } 15.5-17.5 14.5-19.8 15.2-21.0 14.5-19.0 
Salinity 34.28-34.78 | 34.44-34.88 | 34.35-34.81 | 34.50-34.78 
UI TO | 11.0-15.5 10.0-14.5 11.0-15.2 11.0-14.5 
Salinity 34.46-34.60 | 34.40-34.70 | 34.41-34.79 | 34.42-34.62 
I TCO 5.3-10.0: 5.6-11.0 6.4-11.0 
Salinity 34.33-34.64 | 34.34-34.52 | 34.30-34.50 


Table 2.1. Ranges of Temperature and Salinity of Water Masses 
in the northern South China Sea (from Su and Weng, 1994) 


B. | MONSOON EFFECT ON THE CURRENT SYSTEM 
The entire region of the South China Sea is under the influence of the monsoon 
system, and in the absence of major oceanic inflow, the currents undergo a seasonal reversal 
of direction. | 
1. Summer monsoon 
During the summer monsoon season, a warm and weaker southwest wind blows over 
the SCS. The southwest-northeast oriented jet stream is positioned at the central SCS 


(Figure 2.2). Such a typical summer monsoon pattern lasts nearly four months (mid-May to 


mid-September). The summer monsoon forces a wide, uniform northeasterly flowing surface 
current over most of the basin (Figure 2.3). Westward intensification generates a strong 
current along the western side of the basin, off the coasts of Thailand and Vietnam. The 
larger part of the surface flow passes south of Taiwan and merges with the northerly flowing 
Kuroshio; a smaller amount flows north through the Taiwan Strait. Wyrtki (1961) reported 
the existence of a counter current on the eastern side of the basin, with a southwesterly flow 
along the Borneo coast. This flow turns back northerly in the region of Natuna Island, in the 
southern SCS, forming an anticyclonic gyre. Wyrtki hypothesizes that an insufficient supply 
of ae into the basin from the Java Sea produces this gyre. In September, these movements 
decrease and in October the northeast monsoon begins blowing with considerable strength, 
resulting in the complete reversal of the current fields. 

2 Winter Monsoon 

During the winter monsoon season, a cold northeast wind blows over the SCS; as a 
result of the Siberican high pressure system located over the east Asian continent. Radiative 
cooling and persistent cold air advection maintain cold air over the SCS. Tne northeast- 
southwest oriented jet stream is positioned at the central SCS (Figure 2.4). Such a typical 
winter monsoon pattern lasts nearly six months (November to Apmil). The mean northeast 
wind stress over the Western Pacific: deflects the water masses of the North Equatonal 
Current to the south, causing them to enter the SCS through the Luzon Strait. Water from 
the East China Sea moves southerly through the Tarwan Strait and enters the SCS. This ° 


southerly moving water turns southwesterly and flows along the coast of China and Vietnam, 


10 


with a remarkable westward intensification (Wyrtki, 1961). In the center of the SOs: a 
northerly flowing counter current develops over about 10 degrees of latitude (Figure 2.5). 

During the early period of the northeast monsoon, not all of the southerly flowing 
water can flow into the Java Sea, since the southeast monsoon is still blowing weakly over 
this region. This situation produces a northward flowing counter current along the coast of 
Borneo. The water flows as far north as Luzon, where it turns southward and merges into 
the main current, forming a closed cyclonic gyre. 

During December the northwest monsoon is fully developed and the currents are their 
strongest, exceeding 100 cm/s off the coast of Vietnam. The northward flowing counter 
siren in the central SCS is located farther to the east, due to the strength of the main 
current. A weak southward flowing coastal counter current forms along the northwest coast 
of Eee! 
| G; MESOSCALE EDDIES 

Seasonal occurrence of the SCS eddies have been reported by several authors. 
Huang and Wang [1994] mentioned that eddies in the SCS are predominantly cyclonic in 
winter and anticyclonic in summer, with sizes from small (100 km) to medium scale (500 
km). Dale [1956] and Uda and Nakao [1972] reported a cold eddy off the central Vietnam 
coast in summer. Nitani [1970] found a cold eddy located northwest of Luzon in summer. 
Reports from the South China Sea Institute of Oceanology (SCSIO) [1985] indicated that a 
warm-core eddy appeared in summer and winter in the central SCS, but in summer, it was 


closer to Vietnam at the surface. Recently, a cold-core eddy was detected in the central SCS 


I 


between December 29, 1993, and January 5, 1994, from the analysis of TOPEX/ POSEIDON 
data [Soong et al. 1995]. Chu et al. [1997a] and Chu and Chang [1997] identified the 
existence of a central SCS surface warm-core eddy in mid-May from a historical data set: the 
U.S. Navy's Master Observational Oceanographic Data Set ©MOODS). From the composite 
analysis of the U.S. National Centers for Environmental Prediction (NCEP), monthly sea 
surface temperature (SST) fields (1982-1994), Chu et al. [1997b] found that during the 
spring-to-summer monsoon transition (March-May) a warm anomaly (greater than 1.8°C) 
was formed in the central SCS at 112°-119°30'E, 15°-19°30N. 
D. AN AIR-SEA FEEDBACK SCENARIO 

The circulation at intermediate to upper layers of the SCS are mainly influenced by 
monsoons. The prevailing summer (winter) monsoon results in an cyclonic (anti-cyclonic) 
circulation of the Rees layer. The reason may be that the SCS is large, semi-enclosed and 
only slightly affected by the Kuroshio [Huang et al.,1994]. Chu and Chang [1997] 
investigated the historical SST data from the MOODS, and proposed an air-sea feedback 
scenario for the central SCS warm-core eddy formation (Figure 2.6). The anticyclonic surface 
winds prevail in the central SCS from late winter to spring. This anticyclonic wind forcing 
generates downwelling in the central SCS and in turn prevents the cold deep water from being 
advected to the surface. This will promote the formation of the warm-core eddy in the central 
SCS. After the summer monsoon enter an atmospheric surface cyclone occupies the central 
SCS, which entrains the deep cold water into the surface mixed layer. This upwelling effect 


finally destroys the central SCS warm-core eddy. 
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Figure 2.1 Distribution of water masses at the sea surface (a-d) and bottom layer 
(e-h) in the northern South China Sea. a, February; b, May; c, August; 
d, November; e, February; f, May; g, August; h, November; 
see text for acronyms (from Su and Weng, 1994) 
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Figure 2.2 Mean Atmospheric Surface Streamline Analysis 
for summer monsoon (From Cheang, 1980) 
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Figure 2.3 Observational Summer (June) surface Circulation of 


the Southeast Asia Waters (From Wyrtki, 1961) 
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Figure 2.4 Mean Atmospheric Surface Streamline Analysis 
for winter monsoon (From Cheang, 1980) 
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Figure 2.5 Observational Winter (December) Surface Circulation of 
the southeast Asian Waters (From Wyrtki, 1961) 


iv 


Central SCS Surface Atmospheric Anti-cyclone 


Central SCS Downwelling 


Isolating Deep Cold Water 


Promotes Surface Warm-Core Eddy 


Atmospheric Low Pressure at the Central SCS 


Atmospheric Monsoon Onset over SCS 


Downwelling Effect Disappears 


Figure 2.6 An Air- Sea Feedback Scenario for the SCS Warm-Core Eddy Life Cycle 
(From Chu and Chang 1997) 


Il. MASTER OCEANOGRAPHIC OBSERVATION DATA SET (MOODS) 


The MOODS is a compilation of observed ocean data worldwide consisting of (a) 
temperature-only profiles; (b) both temperature and salinity profiles; (c) sound-speed profiles; 
and (d) surface temperature (drifting buoy). These measurements are, in general, irregular 
in time and space. The data are from the early nineteen century to the present and include 
189,059 temperature profiles (Figure 3.1) and 14,004 salinity profiles (Figure 3.2). In this 
study, we analyze the temperature profiles obtained during 1968-84 from a variety of 
instruments. Our study domain includes the area S°N to 25°N and 105°E to 125°E. The 
primary editing procedures included removal of profiles with obviously erroneous location, 
profiles with large spikes, and profiles displaying features that do. not match the characteristics 
of surrounding profiles. In shallow water, this procedure can be partially automated, but it 
also involves subjective interpretation because of the under sampling of MOODS, compared 
to the spatial and temporal variability of the ocean. 

The temporal and spatial distribution of the MOODS data are irregular. The temporal 
distribution (Figure 3.3) indicates that most observations were obtained during the Vietnam 
War. There are nearly 14,800 profiles (maximum number of observations in a year) in 1968 
(Figure 3.4). The minimum number of observations was in 1984 (near 3,00 profiles). There 
figures indicate that the number of observations within a given year exhibits a seasonal 


variation. 
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The main limitation of the MOODS data set is its irregular distribution in time and 
space. Certain periods and areas are over sampled while others lack enough observations to 
gain meaningful insight. Vertical resolution and data quality are Ais highly variable 
depending a great deal on instrument type and sampling expertise. 

We binned seventeen years’ temperature profiles (from 1968 to 1984) into monthly 
data for an individual year which produced 204 temperature data sets. For the salinity 
profiles, we binned into monthly climatology data bins which produced 12 mean monthly 


salinity data sets. 
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Figure 3.1. The Spatial distribution of MOODS temperature profiles 
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Figure 3.2 The Spatial distribution of MOODS salinity profiles 
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Figure 3.4 The temporal distribution of MOODS in 1968 
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IV. ESTABLISHMENT OF GRIDDED DATA 


For each temperature (or salinity) profile, we linearly interpolated the data into 
thirteen standard depths (Table 4.1). After this process, we used the minimum curvature 
spline method [Smith and Wessel, 1990] from the Generic Mapping Tools (GMT) [Wessel 
and Smith, 1992] to create horizontal gridded data. This method interpolates the data 9, 
© at observational points (x, y;) omens grid points with a spline function d(x) such that 


at observational points; 


b(x; 5 ¥;) = 95 | (0) 


The selection of the spline function is based on. the minimum curvature principle, 


which requires continuous second derivatives and minimal total squared curvature, 


C= ff (Poy & d 0) 


Briggs [1974] showed that minimizing C is equivalent to solving the differential equation, 


VV) = DF, S(x-x, + 9-y;) | (2) 


where f , must be chosen such that Equation (1) should be satisfied and where 8 is the 


Kronecker delta function. To solve Equation (2) we need the following boundary conditions 
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ob 


- (3) 
dn? 
and 
0 
—(V*p) = 0 (4) 
On 
along the edges, where o/on indicates a derivative normal to an edge, and 
Z 
TP = 9 (5) 





ox oy 


at — comers. Equations (3),(4),(5) are called free edge conditions; and with these 
conditions, equation (2) has a unique solution with continuous second derivatives called the 

natural bi-cubic spline. The nomenclature a from an analogy with elastic-plate flexure. 
Small displacements of a thin elastic plate of constant flexural rigidity D, subject to a 
vertical normal stress q and constant horizontal forces per unit vertical length of T,,, T,,, T,,, 


approximately satisfy 





Zz Z , 
DV?(V?) - tee 2 pe T,,—] = (6) 
ax? ay oxoy a oy 2 

[Love. 1927]. The minimum-curvature gridding equation (2) is a special case of equation 
(6) when horizontal forces are zero, and the boundary conditions represent zero bending 
moment on the edges (3), zero vertical shear stress on the edges (4), and zero twisting 


moment at the corners (5) [Timoshenko and Woinowsky-Kneger, 1968]. 
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Suppose that T,,= T,, = T and T,, = 0, then Equation (6) becomes 


DV (Vo) - TV’ = q (7) 


When T = 0, Equation (7) is equivalent to Equation (2); but for an arbitrarily large T, the 
solution is dominated by the second term. Here, T has units of force per unit length and the 
T required to adjust the solution scales with D and q; Smith and Wessel [1990] avoided this 
problem by writing 


(1-T)V?(V) - TV) = DF,8(x-x; . y-y;) | (8) 


where 7, is a tension parameter and the / subscript indicates internal tension. When 7, = 0, 
equation (8) reduces to equation (2); and therefore, the minimum-curvature solution is one 
Pidmember case of equation (8). When 7, = 1, the first term in equation (8) vanishes; and 
the solution is harmonic between constraining points. One may prefer to think of this end 
member as representing the steady-state temperature field in a conducting plate with heat 
sources or sinks at the data points. For any 7, in 0< 7,< 1, equation (8) gives a solution with 
continuous curvature, although it does not minimize-equation (1) except when 7, = 0 [Smith 
Pid Wessel, 1990]. 

For the boundary condition, they implemented conditions (4) and (5) but replace 


condition (3) by 





(1 9, 28 wa a0 (9) 


dn? * an 


Za 


where 7, is a tension parameter for the boundary which also varies between 0 and 1. The 
free-edge condition corresponds to 7,;= 0; 7, = 1 forces the solution to flatten at the edge 
(Figure 4.1a,b). This flattening is sometimes desirable, as when gridding potential anomalies 
which should decay toward a regional background field away from the source region. 
Gnidding with tension is an improvement because it adds a degree of freedom and relaxes the 


minimum-curvature constraint. 





Layer Depth (m) Depth (m) 
es 


SS 
ps0 Taso 
allies aaa 





Table 4.1 = standard depth of 13 layers. 
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a) 





Figure 4.1a Cross-sections through surfaces produced with splines 
in tension. The black squares are data constraints. The heavy line 
is the minimum-curvature end member, the thin line is the 
harmonic end member, and the dashed line is an intermediate 
case using some tension (from Smith and Wessel, 1990) 
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Figure 4.1b Cross-sections through continental shelf and slope. The black 
squares represent the intersection between the measured bathymetry (dashed 
line) and 100m isobath contours. These intersections (contour coordinates) were 
then gridded using minimum curvature (heavy solid line) and some tension (thin 
solid line) The minimum-curvature method introduces an extraneous inflection 
point and exceeds the -100m level (The bathymetry in this region is bounded by 
the -200m and -100m level). The surface produced with tension gives a much 
better approximation since it suppresses local maxima and minima between data 
constraints (from Smith and Wessel, 1990) 
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V. COMPOSITE ANALYSIS 


Our goal is to investigate the temporal and spatial variabilities of the SCS thermohaline 
structure. We use Y(X,,),,Z,, T),t,,) to represent temperature and salinity gridded data, where 
(x; ,y) are the honzontal gnds, z, (= 1, 2,..., 13), is the vertical level, 1, 
(=1968,1969,...1984) is the time sequence in years, andz, (€ 1, 2,...,12) is the time 


sequence in month. We define the following temporal average: 


— 1 
W(X VjrZpoF,) = as W(X; jrZp0T bq) (10) 


which 1s the long-term mean value for the month #,, and 


me te _ cM 
W(X;5¥;Z,) = pe W(X iV irZpob,)9 (11) 


which is the ensemble (climatological annual) mean value. Here, at = 17 years. 
A. ANNUAL MEAN 


The annual mean (1968-1984) temperature 7 field over the SCS shows the pattern 


_ of a ee oriented isotherms at the upper layer from the surface to 75 m depth 
(Figure 5.1a). The annual surface mean has a rather weak horizontal temperature gradient, 
decreasing from 28.5°C near the Borneo coast to 25°C near the southeast China coast. A 
stron temperature front is found near the Luzon Strait in the sub-surface layer from 100 m 


to 400 m (Figure 5.1b). The SCS water temperature is quite uniform (e.g., near 15°C at 200 


oa 


m) and much lower than the water east of the Luzon Strait, the West Pacific Ocean Water. 
This might imply that the annual mean Kuroshio intrusion is not significant in the subsurface 
layer. | 

In the central SCS region at 15°N zonal cross-section (Figure 5.2a), the isotherms 


of 7 bend slightly downward from west to east and no frontal structure is found. However, 
in the northern SCS region at 20°N zonal cross-section, the isotherms of 7 show a strong 


frontal structure near 120° to 125°E (Figure 5.2b), which may suggest two distinguished 
water masses, the SCS water mass (both surface and under-surface; «nd the West Pacific 
Ocean water mass. The front between the two water masses strengthens from surface to 400 
m depth. 

The annual mean salinity 5 shows more complexity than temperature. In the upper 
layer (Figure 5.3a), a large amount of fechearen enters the SCS from the Zhujiang (Pearl) 
River in the northwest and from the Me kone River in the southwest. The Kuroshio brings 
the salty water through the Luzon Strait into the north SCS and forms a salty tongue (34 psu) 
stretching into the southeast China coast. As the depth increases, both the Kuroshio intrusion 
and the river run-off effects become less wn At 75 m depth, the Kuroshio intrusion 
effect is still evident, but not the river run-off effect. Below 200 m depth (Figure 5.3b), the 
salinity is more uniform (near 34.6 psu at 200 m) throughout the whole SCS. 

In the central SCS region at 15°N zonal cross-section (Figure 5.4a), the isohalines 


of § are almost parallel to the isobaths in the upper 100 m. Below 100 m depth, the salinity 


a2 


is very uniform (34.4 psu). However, in the northern SCS region at 20°N zonal cross-section 
(Figure 5.4b), the 34.6 psu isohaline clearly shows the Kuroshio intrusion from the Luzon 
Strait into the northern SCS in the layer between 100 and 220m reaching 115°E. 

B. SEASONAL VARIATION 


The long-term monthly mean values relative to the annual mean 


P(X jVjo% ob) = W(x; 59549) s U(x;, 5924) (12) 


represent the composite features of the monthly mean variations of temperature and salinity. 
| Temperature 
a. Surface 
We choose the SST to represent the thermal feature of the surface layer. ° The 
seasonal SST en obtained from the MOODS data is quite consistent with earlier 
investigations based on the National Center for Environmental Predictions (NCEP) data [Chu 
et al., 1997b]. 


During winter (January), T is negative almost everywhere throughout the 


whole SCS (Figure 5.5a). The typical winter (January) pattern contains 
northeast-to-southwest oriented isotherms in the northern SCS (north of 12° N), a warm 
anomaly a > -1.2°C) in the west of Borneo- Palawan islands (WBP), and a cool anomaly 
(T < -2.4°C) near the South Vietnam coast (SVC) in the southern SCS (south of 12°N). The 
spring-to-summer transition (March to May) pattern shows the northward expansion of the 


WEP warm anomaly and the formation of a central SCS (10°-15° N) warm anomaly with ia 


5S 


>0.6°C. The isoline 1°C encloses almost all the SCS in May. 

The summer (July) pattern is characterized by northeast-to-southwest oriented 
isotherms in the northern SCS (north of 15°N) (Figure 5.5b), a WBP cool anomaly (7 < 
0.4°C), and a SVC warm anomaly (> 1.2°C) in the southern SCS. The summer pattern is 
opposite to the winter pattern. The fall-to-winter transition (October) pattern shows the 
northward expansion of a WBP cool anomaly and the formation of a central SCS (10° - 
15°N) cool anomaly. 

b. Sub-Surface Level (75 m) 

We chose 75 m depth to represent ime sub-surface (Figure 5.6). The seasonal 
variation of this layer is characterized by the west Luzon coast (WLC) and SVC cool (warm) 
anomalies and southeast China coast (SCC) warm (cool) anomaly in winter (summer). The 


WLC cool anomaly appears in December, reaching its maximum strength in January with T 


< -3.4°C (Figure 5.6a), lasting two months and disappears in April. The WLC: warm 
anomaly forms in June with T >2°C, reaches T > 2.4°C in July (Figure 5.6b), and lasts for 
four months. 

Figure 5.7 shows the zonal cross-section at 19°N of the T field, representing ° 
the vertical structure of the WLC cool eddy in winter (Figure 5.7a), and the vertical structure 
of the WLC warm eddy in summer (Figure 5.7b). Both WLC winter cool and summer warm 
eddies reach 250 m depth. The winter cool-core (T < -3°C) and the su »mer warm-core (fi > 


2.5°C ) are located between 50 m to 100 m depth. 
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é Intermediate Level (200 m) 

At 200 m depth, the seasonal variation of temperature is featured by a basin- 
wide cool pool ( 7'< -0.8°C), a strong SCC warm pool (T> 0.8°C), and a weak WBP warm 
pool GE > 0.2°C) during the winter (Figure 5.8a) and a basin-wide warm pool Ci > 0.4°C), 
and a weak SVC cool pool (T < 0°C) during the summer (Figure 5.8b). The WLC winter 
cool eddy (T <-1.0°C ) and summer warm eddy Ue >0.8°C ) are still quite evident at that 
depth. The weak WBP warm pool strengthens, expands toward the northeast in spring 
(Figure 5.8a), gradually occupies the whole SCS basin, and turns to the summer pattern in 
July (Figure 5.8b). In autumn, the weak SVC cool pool strengthens and expands toward 
a gradually occupies the whole SCS basin and turns to the winter pattern in January. 
Below 300 m depth, there is no significant seasonal variability. 

2. Salinity 
The seasonal salinity variation § is evident only in the layer above 200 m depth. 
Below this depth, the seasonal variability is very weak. 
a. Surface 


During winter (January), S is positive almost everywhere throughout the 


whole SCS except the SVC and the Gulf of Tonkin (Figure 5.9a). A salty tongue withS >0.4 
psu stretches from the Luzon Strait to the continental shelf along the southeast China coast, 
representing a strong Kuroshio intrusion. Another weak fresh pool withS < -0.2 psu is found 


to the west of Luzon Island. During summer (July), a salty tongue with S around 0.4 psu 
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stretches from the Vietnam coast to 118°E between 10°-17°N, occupying the central SCS 
(Figure 5.9b). Everywhere else is occupied by negative §. A strong fresh tongue occupies 
the continental shelf along the southeast China coast and the Gulf of Tonkin. The negative 
values in the vicinity of the Luzon Strait indicate the cease of the Kuroshio intrusion into the 
SCS. On the other hand, the negative S in the southern SCS shows the influx of fresh water 
from the southern boundary. 

b. Sub-Surface Level (75 m) 

The seasonal variation is featured by a wide salty tongue stretching 
southwestward from the Luzon Strait to the SCS basin in the winter (Figure 5.10a) and a 
fresh tongue stretching northward from the southwest boundary near the Mekong River 
mouth to the SCS basin in the summer (Figure 5.10b). The former represents the intrusion 
of the salty Kuroshio _— (with maximum iS _ 0.6 psu) through the Luzon Strait, and the 
latter indicates the low salinity water (with minimum §<-0.4 psu) extending to the SCS. 

Figure 5.11 shows the zonal cross-section at 19°N of the S field, representing 
the vertical structure of the WLC high salinity anomaly (S > 0.4 psu) in winter (Figure 5.11a) 
and low salinity anomaly (S < -0.2 psu) in summer (Figure 5.11b). Both WLC high Jow 
salinity anomalies are located between 50 m to 100 m depth. 

C. Intermediate Level (200 m) 
At 200 m depth, the seasonal variation of §: shows a similar pattern (Figure 
> 12) to that at 75 m depth. In the winter (Figure 5.12a), positive § with a maximum value 


of 0.3 psu occupies the most area of the SCS, except the southwest corner near the Mekong 
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River, where § is negative (-0.4 psu). In the summer (Figure 5.12b), a narrow fresh tongue 


with a minimum value of -0.4 psu stretches northeastward from the southwest corner near the 


Mekong River to 20°N (Figure 5.12b). 
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Figure 5.la The annual mean temperature field (7) 
from surface to 100 meters 
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Figure 5.1b The annual mean temperature field (7) 
. from 125 to 400 meters 
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Figure 5.2a The zonal cross-section view of annual mean temperature 7 
from Latitude 14°N to Latitude 19°N 
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Figure 5.2b The zonal cross-section view of annual mean temperature 7 
from Latitude 20°N to Latitude 25°N 
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Figure 5.3a The annual mean salinity field (S$) 
from surface to 100 meters 
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Figure 5.3b The annual mean salinity field (5) 
from 125 to 400 meters 
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Figure 5.4a The zonal cross-section view of annual mean salinity (S$) 
from Latitude 14°N to Latitude 19°N 
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Figure 5.4b The zonal cross-section view of annual mean salinity ($) 
from Latitude 20°N to Latitude 25°N 
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Figure 5.5a Monthly mean temperature anomalies (7) at surface 
from January to June 
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Figure 5.5b Monthly mean temperature anomalies (7) at surface 
from July to December 
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Figure 5.6a Monthly mean temperature anomalies (7) at 75 meters 
from January to June 
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Figure 5.6b Monthly mean temperature anomalies (7) at 75 meters 
from July to December 
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Figure 5.7a Zonal cross-section view of monthly mean temperature 
anomalies (7) at 19°N from January to June 
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Figure 5.7b Zonal cross-section view of monthly mean temperature 
anomalies (7°) at 19°N from July to December 
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Figure 5.8a Monthly mean temperature anomalies (7) at 200 meters 
from January to June 
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Figure 5.8b Monthly mean temperature anomalies (T ) at 200 meters 
from July to December 
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Figure 5.9a Monthly mean salinity anomalies (S ) at surface 
from January to June 
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Figure 5.9b Monthly mean salinity anomalies (5) at surface 
from July to December 
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Figure 5.10a Monthly mean salinity anomalies (5) at 75 meters 
from January to June 
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Figure 5.10b Monthly mean salinity anomalies (S) at 75 meters 
from July to December 
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Figure 5.11a Zonal cross-section view of monthly mean salinity 
anomalies (5) at 19°N from January to June 
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Figure 5.11b Zonal cross-section view of monthly mean salinity 
anomalies ($) at 19°N from July to December 
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Figure 5.12a Monthly mean salinity anomalies (5) at 200 meters 
from January to June 
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Figure 5.12b Monthly mean salinity anomalies (5) at 200 meters 
from July to December 
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VI. EMPIRICAL ORTHOGONAL FUNCTION (EOF) ANALYSIS 


The empirical orthogonal function (EOF) analysis is a useful tool to describe the main 
mode of temporal and spatial variabilities. It is usually applied to the anomaly data relative 
to the climatological mean. 

A: THERMAL VARIABILITY 


The temperature anomaly field relative to the climatological annual mean field 


NEE ecstnt,) = Neesy scsi TCV acer) : (13) 


is re-arranged into a NP eae T(r, ; i) ,n=1,2,...,N; and p=1, 2,..., P. Here P 
= 204 1s the total number of time points used for computing the covariance matrix, that is, 17 
years of monthly data sets; N = $733, scccosiiill to the number of Onids\(ta ls 2 ly 
= 1, 2,...,21; k= 1, 2,..., 13). The non-seasonal thermal variability is investigated using the 
Empirical Orthogonal Function (EOF) analysis. This method separates the data sets into 
eigenmodes. Generally speaking, each mode has an associated variance, a dimensional spatial 
‘ pattern, and a non-dimensional time series. For example, Lorenz (1965) used the eigenvectors 
of the covariance matrix associated with a sample data field to calculate the time series of 
coefficients or predictors for the statistical weather prediction. Chu et al. (1997a, b) used this 
method to investigate the non-seasonal variabilities of SST and surface wind stress. From 


our data 7 (ie, t,) , a 5733x5733 spatial covariance matrix is calculated by 
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a Ral 
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DY Trt) Tryst)» N=5733, P=204 (14) 
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ao 

R 
Dea ose 

where 7 and m (1,2,...,5733) denote the grid locations, The diagonal elements of the 

covariance matnix R,,, =(7=1,2,...,N) are the variance at location 7,, The off-diagonal elements 

are the covariance with spatial ae equal to the difference between the row and column 

indices. This symmetric matrix has its 5733 real eigenvalues A, , and eigenvectors ©, (r;) 


such that 


y 1200 (a) OM eee = The al ce) 
i 

The eigenvectors ®,, ®,, ..., ®y are called Empirical Orthogonal Functions with the unit of 
°C. Each ®, is a 5733-point (212113) distribution of i variation temperature 
anomaly pattern. The eigenvalues, A,(a=1, 2, ..., N), are all positive and the summation of 
them, ~A,, equal the total variance. Therefore, A, is considered as portion of total variance 
"explained" by the EOF ®,. It is convenient to label the eigenfunctions, ®,, so that the 


eigenvalues are in descending order, i.e, A, > A, > A, >... 


The data matrix 7 (Tire t,) is thus written approximately by 


Pee, opt > PCC potg Py (X;2¥ 024) (16) 
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where PC,(t, ¢,) is the principal component representing the temporal variation of the 
associated spatial pattern described by the EOF mode, ,(, y,,z,). 
B. PRINCIPAL EOF MODES 

The first six leading EOF modes are able to account for 64.1% of the total variance 
(Table 6.1). The first two leading modes, EOF1 and EOF2, explain 44.4% of the total 
variance. Each EOF mode is three-dimensional and normalized so that its total spatial 
variance is equal to unity. We may easily see the vertically coherent structure of each mode. 

ae The EOF1 Mode 

The EOF! mode (@,) shows frontal and multi-eddy structures. The frontal structure 
is evident in the northern SCS. Only positive values occur in the upper layer (above 75 m 
depth) and both positive and negative values appear at the other depths (Figure 6.1a,b). 
However, we cannot use these positive/negative values of @, to identify the warm or cool 
anomaly of the individual eddy. This is because the temperature anomaly, due to the EOF 1 
pattern, is the product of PC, and @,. 

a. Western Bena Current and Frontal Structures 
From the surface to 30 m depth, a bi-frontal structure (northeast-to-southwest 

oriented isolines) is the major thermal feature in the northern SCS (north of 12°N): a 
strong coastal front along the China coast, and a relatively weak and wide front across the 
SCS basin from the Vietnam coast to the Luzon island between 12°-17°N. The two fronts 
change their strengths synchronously, since the values of ®, are positive for the two fronts. 


The bi-frontal structure weakens with the depth and disappears at 50 m depth. The front 
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reoccurs from 75 m depth and continues to be seen down to 300 m depth. Below 75 m 
depth in the northern SCS, @, shows a northeast-to-southwest oriented single front structure 
aligning to the China coast. The values of 9, across the front change signs from negative 
near the China coast to positive in the middle of the northern SCS, indicating the out-of- 
phase variability of the front near the China coast and the deep basin. 
Furthermore, we may notice that the values of @, change signs from the top 
layer (0-30 m) to the lower layer (100-300 m) near the western boundary. This implies a 
strong baroclinicity of the western boundary currents with the opposite directions at the top 
layer and the lower layer. 
b. Southern SCS Multi-Eddy Structure 
In the southern SCS (south of 12° N) from the surface to 30 m depth, the 
EOF 1 — a dipole structure (dual eddies) with a high center (0.3°C) near the mouth of 
the Mekong River, and a low center (< 0.14°C) in the west of the Borneo-Palawan islands 
(WBP). This dipole structure is still quite evident from 30 m to 75 m depths. Below 75 m 
depth, the western eddy, near the Mekong River mouth, is much stronger that the eastern 
- eddy. 
| At 30 m depth, an eddy (©0.4°C) appears in the northern SCS. This eddy 
slants towards the east (the Luzon Strait) with depth and becomes very strong in the layers 


between 75 m and 150 m. Below 150 m depth, this eddy reduces its strength. 
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C. Kuroshio Intrusion 
In the northern SCS, the isolines of @, intersect the cross-section of the 
Luzon Strait with the patterns showing Kuroshio intrusion/SCS outflow oon the Luzon 
Strait. This pattern indicates a weak Kuroshio intrusion/SCS outflow in the upper layer 
above 50 m depth, becomeing stronger below 50 m depth and reaching maximum strength 
in the layer between 100m and 150 m depths. 
Ds The Second EOF Mode 
Similar to the first EOF mode, the second EOF mode (@,) has only positive values 
in the upper layer from the surface to 50 m depth and has both positive and negative values 
appearing on the other depths. The negative @, area is limited in the southeast corner of the 
SCS with small absolute values (Figure 6.2a,b). This means that the overall SCS basin 
thermal variability is either warming or cooling. We may call this mode the 
warming/cooling mode. In addition, the isolines of ®, at all the depths near the Luzon Strait 
are almost parallel to the cross-section of the strait. This indicates that the second EOF 
(warming/cooling) mode is not caused by the Kuroshio intrusion/SCS outflow. It might be 
related to the atmospheric forcing. Furthermore, the second EOF mode also shows eeratal 
front and multi-eddy structures. 
a. Coastal Front Structure 
In the north SCS, a coastal front structure can be identified on the continental 
shelf near the China coast (Figure 6.2a). As the depth increases, this front structure extends 


to the southern SCS and occupies the western part of the SCS from 100m to 300 m depths 
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(Figure 6.2b). This strong gradient of ®, may be caused by the western boundary current. 
b. Northern SCS Multi-Eddy Structure 
The EOF2 clearly shows a multi-eddy structure in the northern SCS sub- 
surface layer between 50 m and 125 m. At 50 m depth, a strong north-south dipole structure 
(dual eddies) with a high center (0.2°C) in the west of the Luzon Strait (18°-22°N, 115°- 
119° FE), and alow center (< 0.04°C) in the central SCS (12°-15°N, 110°-116°E) . This 
dipole structure strengthens at 75 m depth and weakens as depth increases from 75 m. 


Below 125 m depth, the northern SCS multi-eddy structure disappears. 





iano 





Table 6.1 cease of the First Six Leading EOFs 


C. TEMPORAL VARIABILITIES 
Each principal component can be treated as the projection of the temperature 


anomaly through a ‘filter’ of corresponding EOF mode (Chu et al., 1997). 
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i: First Principal Component 


The first principal component, PC,(t, , ¢,) , 18 the time series (1968-1984) of 


the EOF 1 amplitudes and shows basically a seasonal variation. During the winter (summer) 


monsoon season, PC, (1, ) has negative (positive) values (Figure 6.3), which indicate 


that the temperature anomaly, due to this mode, has the opposite (same) sign with respect 
to EOF1. Since the first EOF mode @, is always positive throughout the upper layer (0 - 


50 m) of the whole SCS (Figure 6.1), positive (negative) values of PONG ) during 


summer (winter) monsoon season correspond to warm (cool) anomalies. This shows the 
regular seasonal variation. 
An interesting feature of the first EOF mode @, is its negative value near the western 


boundary in the sub-surface layer (75 - 300 m). Since PC, (1, ) has negative (positive) 


values during winter (Summer) monsoon season, the thermal variability of the western 
boundary is out-of-phase with respect to the seasonal variability, that 1s, warm (cool) 
anomaly in the winter (summer) monsoon season. The warm anomaly in the winter could 
be related to the Kuroshio intrusion. | 


Accordingly, the seasonal variation represented by PC, ( t, ) is not a pure 


sinusoidal. The summer had a larger variability in the 1980's than in the early 1970's. 


However, the winter had a larger variability during 1968-1972 than the rest of the penods. 
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2. Second Principal Component 


The second principal component, Be: (t, ) _, is the time series (1968-1984) of the 


second EOF (warming/cooling) mode amplitudes and shows interannual variabilities with 
periods of two to five years (Figure 6.3). The positive (negative) values of 


PC, (t,) indicate the SCS warming (cooling). The maximum value of 9, is 0.3°C near 


the China coast at the surface and a high center with isoline of 0.2°C in the west of Luzon 
Strait at depths between 50 m arc 125 m (Figure 6.2). 


The values of PC, ( t, ) vary between 10 and -10. Its maximum value (10) 


appeared in August 1969 and September 1976, which indicates the strong SCS warm 
anomaly with the maximum warming of 3°C [0.3°Cx10] at the surface near the China coast, 
and the occurrence of a sub-surface (50-125 m) warm eddy with the anomaly of 2°C 
[0.2°Cx10] in the west of the Luzon Strait. Its minimum value (-10) appeared in October 
1974, which indicates the strong SCS cool anomaly with the maximum cooling of -3°C 
[0.3°Cx(-10)] at the surface near the China coast, and the occurrence of a sub-surface (50- 
125 m) cool eddy with the anomaly of -2°C [0.2°Cx(-10)] in the west of the Luzon Strait, 

We found the following criterions in evident warm/cool anomaly occurrence with the 


maximum temperature anomaly exceeding 1.5°C. When PC,(t, ) >5, the SCS has an 


evident warming (maximum 1.5°C) and a warm eddy occurs at the west of the Luzon Strait 


in the sub-surface with the warm anomaly of 1°C. Onthe other hand, PC,( t, y) =>) te 
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SCS has an evident cooling (maximum above 1.5°C) and a cool eddy occurs at the west of 
Luzon Strait in the sub-surface with the cool anomaly of -1°C. 

The strong anomaly periods (maximum above 3°C) and evident periods (maximum 
above 1.5°C) are listed in Table 6.2. We may find that the SCS warm anomaly generally 


appeared during the period of 1968-1972. After 1972, the cool anomaly occured more often, 


except in September of 1976. 


Anomaly Maximum Anomaly Maximum Anomaly 
Period Above 3°C Above 1.5°C 


SCS Warming Aug. 1969 Sep. 1968, 
Sep. 1976 Jan. Feb. May Jul. Aug. 1969 

Nov. 1969 

, Jun. 1970 : 

a) | Oct. Dec. 1971 

May Jul. Sep. Nov. 1972 

Sep. 1976 
Jul. 1977 





SCS Cooling Oct. 1974, Oct. 1974 
Jan. 1976 
Feb.- Mar. 1977 
Feb. 1980 
Apr. 1983 


Table 6.2 SCS warm/cool anomaly periods 
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Figure 6.l1a EOF 1 Temperature anomalies (Unit 0.01°C) 
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Figure 6.1b EOF 1 Temperature anomalies (Unit 0.01°C) 
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Figure 6.2a EOF 2 Temperature anomalies (Unit 0.01°C) 
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Figure 6.2b EOF 2 Temperature anomalies (Unit 0.01°C) 
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Figure 6.3 The first six principal components (PC) for 1968-19384 
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VIL. THE SOUTH CHINA SEA CIRCULATION 


A. GEOSTROPHIC VELOCITY INVERTED FROM DENSITY FIELD 
The thermal wind relation is often employed to calculate geostrophic velocity from 
hydrographic data 


oo 
u=u,+ 2 f 
fPo 


> 


dz’ (17) 
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where (u, v) and (u,, v,) are the geostrophic velocity at any depth z and at a reference depth 
Z,» is the in-situ water density, p, is the characteristic value of the density, f= 2QsinA 
is the Coriolis parameter, where Q is the earth's rotation rate and A is the latitude. The 
hydrographic data determine the geostrophic vertical shear. The reference velocity (u,, v.) 
can be determined by several inverse techniques, such as the B-spiral method [Stommel and 
Schott, 1977; Schott and Stommel, 1978], the box inverse method [Wunsch, 1978], and the 
Bernoulli method [Killworth, 1986]. Wunsch [1994] proposed a generalized B-spiral and 
geostrophic balance inverse method based on a nonlinear optimization on Needler's P 
equation [Needler, 1967]. Davis [1978] pointed out that the B-spiral method and the box 


inverse method, regardless of appearances, are based on the same dynamical assumptions, 


and the existence of the B-spiral thus becomes a required condition for these inverse 


dd 


methods. 

Non-geostrophic inverse methods were also developed. Bogden et al. [1993] 
estimated the tume-averaged velocity field in the North Atlantic from observations of density, 
wind stress and bottom topography based on the assumption of minimum mixing of potential 
density at mid-depth. Tziperman et al. [1992] applied the adjoint method to investigate the 
steady-state oceanic general circulation using a realistic numerical oceanic general circulation 
model and hydrographic data. 

Recently, Chu [1995] and Chu et al. [1998a] proposed a simple P-vector inverse 
method to obtain the absolute velocity from hydrographic data. A velocity vector 1s 
represented into the product of a unit vector (P-vector) and a scalar, | 

V = r(x, y, z) P. (19) 
The inversion of V is fulfilled by two steps: (1) determination of the unit vector P, and 
(2) determination of the scalar r(x,y,z). Since the geostrophic flow is along both the 


isopycnal and the iso-potential vorticity surfaces, the velocity V is parallel to Vq x Vp. Here, 


o is the potential density, and 
q=fop/ o | (20) 


is the potential vorticity. If the p surface is not parallel to the q surface, the unit vector P cz. 


be defined by [Chu, 1995] 


_ Vp x Vq 


= (21) 
|\Vp x Vq| 
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Applying the thermal wind relation at two different depths z, and z,,, as shown in 


Figure 7.1, a set of algebraic equations for determining the parameter r is obtained 


k) p(k) (mn) _ | 
r ( Ms aoe - PoP. = Au, (22) 
(k) p©®) _ Lim p(™)_ 
r ie r Lie = Ay, ~ 


which are two linear algebraic equations for r and r™ [r © = r(x,y,z, )]. Here 


g 7 a dp 
Au 4 Av — —— = 9 _ d / 
(Atine Aim) = i (Gy ay (23) 


where 9 is the in-situ density. 

Chu et al. [1998b] shows the advantages of using the P-vector method, namely, the 
easy identification of the two necessary conditions for the B-spiral inverse method: (1) 
intersection between the p- and q-surfaces {existence of the vector P), and (2) the B turning 
(turning of the vector P). 

For water columns satisfying the two necessary conditions, equation (22) may be 
used to compute r™ for the ee z,. There are (N-1) sets (m = 1, 2, k-1, k+1, ... N) of 
equations (22) for calculating r™. Here, N is the total number of vertical levels of the wafer 
column. All of the (N-1) sets of equations are compatible under the thermal wind constraint 
and should provide the same solution. However, due to errors in measurements 
(instrumentation errors) and computations (truncation errors), the parameters r“ may vary 
with m. A least square error algorithm [Chu et al., 1998b] is used to minimize the error. 


For further details and validation of the algorithm presented, see Chu et al. [1998a]. 
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B. SEASONAL VARIATION OF SCS CIRCULATION 

1 General Description 

The most obvious features of the inverted summer and winter SCS circulation are 
flow separation, basin gyres, and mesoscale eddies. Hinted in Wyrtki's (1961) depiction 
but more explicit in our computation, are the upper level (10 m depth) west-to-east 
cross-basin currents, which are almost parallel to 16°N latitude from February to July. The 
inverted monthly SCS circulation agrees quite well with earlier observational study (Wyrtki, 
1961). As the commonly used inverse methods, the P-vector inverted western boundary 
currents are somewhat small. This could be caused by the neglect of relative vorticity in the 
potential vorticity calculation [see (21)]. 

oe Surface Circulation 

During the winter monsoon period (November to March), the winter Asian high 
pressure system brings strong winds from the northeast (Figure 2.4) on the SCS surface 
circulation pattern is cyclonic (Figure 2.5). At the surface inflow from the Luzon Strait (the 
Kuroshio intrusion) turns southwest along the Asian continental shelf, then south along the 
coast of Vietnam (Figure 7.2). In November, this current splits into two currents at 12°N: 
a southward alongshore current and an eastward off-shore current. The southward branch 
turns east at 6°N, then turns northeast and joins the eastern branch at 115°E. An evident 
multi-eddy structure is obtained with two cyclonic eddies in the north SCS (north of 12°N) 
and cyclonic and anticyclonic eddies in the southern SCS (south of 12°N). The flow pattern 


does not change vary much from November to December, except the cyclonic eddy near the 
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Luzon Strait strengths. In January, the west cyclonic eddy in the northern SCS expands 
toward the south. In February, the west cyclonic eddy in the northern SCS disappears and 
instead a cross-basin flow appears along 16°N latitude and joins the cyclonic eddy near the 
Luzon Strait. The southern SCS is occupied by an anticyclonic eddy. This anticyclonic eddy 
expands to north and forms the only evident eddy in the central SCS in March (Figure 7.2a). 
This eddy is associated with the central SCS warm pool in spring (Appendix A) and expands 
toward the south in April. May is the month of the summer monsoon transition. The 
anticyclonic eddy reduces its strength and becomes unidentifiable. 

During the summer monsoon period (Mid-May to August), winds blow from the 
southwest (Figure 2.2) and the SCS surface circulation generally follows suit with 
anticyclonicity in the southern basin (Figure 2.3). Inflow is through the southern boundary 
and outflow is through the northern Taiwan Strait and eastern Luzon Strait. The inverted 
July means the general circulation pattern has the following features. Velocities reach 20 
cm/s near the Vietnam Coast and splits into two currents at 12°N: the coastal current and 
off-shore current. The off-shore current further bifurcates and partially leaves the coast; the 
bifurcation point ts at 110°E. The northern branch moves northeast to 113°E and thet turns 
zonally between 15°-18°N. The southern branch moves zonally until reaching the Palawan 
Island, then turns north and joins with the north branch at 16°N. A cyclonic eddy appears in 


July near the south Vietnam coast and strengthens in August. 
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3 Sub-Surface (75 m) 

The sub-surface (75 m) shows a strong seasonal variation (Figure 7.3). Warm and 
salty Kuroshio water enters the SCS through the Luzon Strait in October-March, the 
transition to and during the winter monsoon season. The winter circulation pattern is the 
basin-wide cyclonic gyre (November - January). A cyclonic mesoscale eddy splits from the 
gyre near the Luzon Strait in January and keeps its strength in February. On the other hand, 
the basin-wide cyclonic gyre weakens and shrinks toward the southwest in February and 
disappears in March. 

During the summer monsoon period (June-September), the Kuroshio intrusion 
through the Luzon Strait ceases. A weak anticyclonic eddy occurs in June near the South 
Vietnam coast and strengthens in July. A cross-basin current establishes in August from the 
north branch of that eddy to the Luzon Strait and forms a weak cyclonic gyre in the 
northwest and an anticyclonic gyre in the southeast. During the winter monsoon transition 
period (October), the southeast anticyclonic gyre weakens and the northeast cyclonic gyre 
gains strength and eventually the cyclonic gyre occupies the whole SCS basin, the beginning 
of the winter circulation pattern. | | | 

4, Intermediate Level (400 m) 

As the water depth increases, the inverted velocity generally reduces with depth to 
200 m. The velocity at 200 m depth is quite small (see Appendix D). However, the velocity 
increases with depth below 200 m. Figure 7.4 shows the monthly variation of the SCS 


circulation at 400 m. At that depth, the inverted flow flows out of the SCS through ‘ie 
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Luzon Strait. The winter pattern is characterized by three anticyclonic eddies (December- 
February), with two eddies in the northern SCS, and one in the southern SCS. The 
anticyclonic eddy near the Luzon Strait is evident from December to February. An east-to- 
west cross-basin current establishes along 16°N in February from the southern branch of 
the Luzon eddy. This cross-basin current sustains from February to April and separates the 
SCS basin into two parts. The northern part is occupied by an anticyclonic eddy and the 
southern part by a cyclonic eddy. The summer circulation pattern is featured by double 
anticyclonic gyres. The southern gyre is quite evident in June and July. The tangential 
velocity is around 4 cm/s. 

= Kuroshio Intrusion Through the Luzon Strait 

The SCS joins the Pacific Ocean through the Luzon Strait. The seasonal variations 
of the intrusion of the Kuroshio Water into the SCS through the Luzon Strait have been 
investigated in earlier studies [Fan and Yu, 1981; Shaw, .1989; Shaw, 1991]. As pointed out 
by Shaw [1991], the intrusion process is important not only to the transport of water masses 
into the SCS, but also to the shelf circulation off the southeast coast of China. Shaw [1989, 
1991] used the discriminant analysis method to classify the water mass T, S characteristics 
at 150, 200, and 250 meters, and found that water characteristics of the Philippine Sea 
(Kuroshio) was identified along the continental margin south of China from October to 
January. The presence of this water indicated an intrusion current from the Philippine Sea 
into the SCS. 


To verify Shaw’s result, we computed the east-west geostrophic velocity along 
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120.5°E longitude using the P-vector method. The temperature and salinity used for the 
computation are the gridded data from the surface to 400 m depth along 120° and 121°E 
(Figure 7.5). The monthly zonal velocity in 120.5°E cross-section 1s given in Figure 7.6. 
The negative values indicate westward flow through the Luzon Strait (Kuroshio intrusion). 
The positive values show the outflow from the SCS to the Philippine.Sea. Flow through the 
Luzon Strait has a strong seasonal variation, as well as a vertical shear. Evident Kuroshio 
intrusion (negative values) occurs during the winter monsoon season (November - March) 
in the upper 200 m. Below 200 m, the velocity is quite small, and most values are positive. 
We also computed the upper layer (0-400 m) volume transport through the Luzon Strait 


(Figure 7.7). The evident Kuroshio intrusion appears from November to March. 
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Figure 7.2b Surface Circulation from J uly to December 
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Figure 7.5 Temperature and Salinity data from surface to 400 meter used 
for the computation of east-west geostrophic velocity along 120.5°E 
(dash lines: T and 8; solid line: geostrophic velocity) 
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Figure 7.62 The monthly zonal velocity in 120.5E from January to June 
( negative value indicate westward flow; Kuroshio intrusion) 
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Figure 7.6b The monthly zonal velocity in 120.5E from July to December 
( negative value indicate westward flow; Kuroshio intrusion) 
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Figure 7.7 The Mass Transport (Sv) at the Luzon strait 
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VOI. CONCLUSIONS 


The goals of this study were to investigate the thermohaline structure and the 
circulation patterns of the South China Sea using the historical MOODS hydrographic data. 
We applied the minimum curvature with spine method to produce thirteen level gridded data 
from the MOODS. Although the MOODS contain more than 180,000 temperature profiles, 
the deeper layer (below 400 m depth) is poorly sampled. Futhermore, the MOODS contain 
only 14,000 salinity profiles, so we were only able to establish climatological annual and 
monthly mean salinity fields at.thirteen levels above 400 m depth. After the establishment 
of the gridded temperature and salinity data, we used the composite and EOF analysis to get 
the seasonal and interannual variabilities of the SCS thermohaline structure, and the P-vector 
inverse method to obtain the dhrco:2itiaalell absolute geostrophic velocity fields. The 
major results from this research are summarized as follows. 

A. CLIMATOLOGICAL MEAN THERMOHALINE FEATURES 


The annual mean (1968-1984) temperature 7 field over the SCS shows the pattern 


of northeast-southwest oriented isotherms at the upper layer from the surface to 20 m depth 
with a weak gradient from 28.5°C near the Borneo coast to 25°C near the southeast China 
coast. A strong temperature front is found near the Luzon Strait in the sub-surface layer 
from 100m to 400m. The SCS water temperature is quite uniform (e.g., near 13°C at 250 
m) and much lower than the water east of the Luzon Strait, the West Pacific Ocean Water. 


The annual mean salinity § shows that in the upper layer, a large amounts of freshwater 


oF 


enter the SCS from the Zhujiang (Pearl) River in the northwest and from the Mekong River 
in the southwest. The Kuroshio brings the high salinity water through the Luzon Strait into 
the north SCS and forms a salty tongue (34 psu) stretching into the southeast China coast. 
B. SEASONAL THERMOHALINE VARIABILITY 

Seasonal thermohaline variability is associated with fronts and multi-eddies. The bi- 
frontal structure is the major thermal feature in the northern SCS (north of 12°N): a strong 
coastal front along the China coast, and a relatively weak and wide front across the SCS 
basin from the Vietnam coast to the Luzon island between 12°-17°N. The two fronts change 
their strengths synchronously; they weaken with the depth and disappear at 50 m depth. The 
front reoccurs from 75 m depth and is evident as deep as 300 m depth. At 30 m depth, an 
eddy (>0.4°C) appears in the northern SCS and slants toward the east (the Luzon Strait) 
with depth and wreceataen very strong in the layer between 75 m and 150 m. Below 150 m 
depth, this eddy reduces its strength. In the southern SCS (south of 12°.N) from the surface 
to 30 m depth, a dipole structure (dual eddies) was found with a warm (cool) eddy near the 
Mekong River mouth, and a cool (warm) eddy in the west of the Borneo- Palawan islands 
during the summer (winter) monsoon season. Bélow 75 m depth, the western eddy near the 


Mekong River mouth is much stronger that the eastern eddy. 
The seasonal salinity variability was also identified, During winter (January), a salty 


tongue 1.2 psu anomaly stretches from the Luzon Strait to the continental shelf along the 


southeast China coast, representing a strong Kuroshio intrusion. A weak fresh anomaly (- 
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0.2 psu) is found in the west of the Luzon island. During the summer (July), a salty tongue 
with 0.4 psu stretches from the Vietnam coast to 118°E between 10°-17°N, occupying the 
central SCS while the rest is occupied by negative anomaly. A strong fresh tongue occupies 
the continental shelf along the southeast China coast and the Gulf of Tonkin. 
G; KUROSHIO INTRUSION 

Both water mass analysis and the inverted geostrophic velocity at the Luzon Strait 
confirm the Kuroshio intrusion ane the winter monsoon season. Such an intrusion has a 
vertical shear. Evident Kuroshio intrusion (negative values) occurs during the winter 
monsoon season (November - March) in the upper 200 m while below 200 m the velocity 
is quite small and outflows from the SCS into the Philippine Sea. The computed upper layer 
(0-400 m) volume transport through the Luzon Strait has a strong seasonal variation: the 
maximum Kuroshio intrusion into the SCS with a value of -2.2 Sv occurmng in March, and 
the maximum SCS outflow with a value of 1.5 Sv in September. 
D. INTERANNUAL SCS WARMING/COOLING PHASES 

The interannual SCS warming/cooling phases were identified by the second EOF 

mode with periods of two to five aes The strongest SCS Bt occurred in August 
1969 and September 1976, with a warming of 3°C at the surface near the China coast and 
an occurrence of a sub-surface (50-125 m) warm eddy (2°C) in the west of Luzon Strait. 
The strongest SCS cooling appeared in October 1974 with a cooling of 3°C at the surface 
near the China coast, and an occurrence of a sub-surface (50-125 m) cool eddy (-2°C) in the 


west of the Luzon Strait. 


92 


The warming events appeared between the late 1960's and first half of the 1970's, and the 
cooling events occurred in the second half of the 1970's and early 1980's. 
E. MULTI-EDDY STRUCTURE 

The T, S and the inverted velocity fields show the multi-eddy structure: two cyclonic 
eddies in the north SCS (north of 12°N) and cyclonic and anticyclonic eddies in the southern 
SCS (south of 12°N). The anticyclonic warm-core eddy expands to the north and forms 
the only evident eddy in the central SCS in March, and expands toward the south in April. 
May is the month of the summer monsoon transition. The anticyclonic warm-core eddy 
reduces its strength and becomes unidentifiable. During the summer monsoon period (mid- 
May to August), winds blow from the southwest and the SCS surface circulation generally 
follows suit with anticyclonicity in the southern basin. Inflow is through the southern 
boundary and outflow is through the arene Taiwan Strait and eastern Luzon Strait. The 
inverted July mean general circulation pattern has the following features. Velocities reach 
20 cm/s near the Vietnam Coast and splits into two currents at 12°N: the coastal current and 


off-shore current. 
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APPENDIX A. MONTHLY TEMPERATURE ANOMALY 
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